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Abstract 

Putative organisms on the Martian surface would be exposed to 
potentially high doses of ionizing radiation during strong solar X-ray 
flares. We extrapolate the observed flare frequency-energy release 
scaling relation to total X-ray energies much larger than seen so far 
for the sun, an assumption supported by observations of flares on other 
solar- and subsolar-mass main sequence stars. Flare spectra are taken 
as power laws, with the logarithmic slope a parameter based on the 
observed statistics of the most energetic hard X-ray flare spectra. We 
calculate the surficial reprocessed spectra using a Monte Carlo code 
we developed for the transport of X-rays and gamma rays, including 
photoabsorption and detailed Compton scattering. Biological doses 
from indirect genome damage are calculated for each parameterized 
flare spectrum by integration over the X-ray opacity of water. The re- 
sulting doses depend sensitively on spectral slope, which varies greatly 
and unsystematically for solar flares. Using the roughly uniform ob- 
served distribution of spectral slopes, we estimate the mean waiting 

*To be published in Planetary and Space Science. 
^ Corresponding author 

■^E-mail: {dss ,parrot}@astro . as .utexas . edu 



1 



Smith 



2 



time for solar flares producing a given biological dose of ionizing ra- 
diation on Mars and compare with lethal dose data for a wide range of 
terrestrial organisms. These timescales range from decades for signif- 
icant human health risk to 0.5 Myr for D. radiodurans lethality. Such 
doses require total flare energies of 10 33 -10 38 erg, the lower range of 
which has been observed for other stars. Flares are intermittent bursts, 
so acute lethality will only occur on the sunward hemisphere during a 
sufficiently energetic flare, unlike low-dose-rate, extended damage by 
cosmic rays. We estimate the soil and CO2 ice columns required to 
provide l/e shielding as 4-9 g cm" , depending on flare mean energy 
and atmospheric column density. Topographic altitude variations give 
a factor of two variation in dose for a given flare. Life in ice layers that 
may exist ~ 100 g cm~ 2 below the surface would be well protected. 
Finally, we point out that designing spacesuits to sufficiently block 
this radiation on Mars missions may be difficult, given the conflict be- 
tween solutions for lightweight protection from energetic particles and 
those from X-rays. 

1 Introduction 

Habitability for all types of life — from human to microbial — is a central problem 
for the next decades of Mars exploration. In particular, knowledge of surficial 
radiation doses from all possible sources is essential. Cordoba-Jabonero et al. 
(2003) have calculated in detail the transfer of ultraviolet (UV) radiation in the 
Martian atmosphere and folded the surface flux with a DNA action spectrum to ob- 
tain UV dose rates (see also Cockell et al., 2000, Mancinelli and Klovstad, 2000). 
Schuerger et al. (2003) treated the same problem experimentally, with the survival 
of spores of Bacillus subtilis in a Mars simulation chamber. Molina-Cuberos et al. 
(2001), Pavlov et al. (2002), and De Angelis et al. (2004) have estimated the fluxes 
of Galactic cosmic-ray particles on the present and ancient Martian surface and 
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subsurface. Pavlov et al. find that the maximum dose rate due to Galactic cosmic 
rays is about 0.2 Gray yr _1 , using dosimetry units in which 1 Gray (Gy) = 100 rad 
and 1 rad = 100 erg g _1 absorbed. This dose rate occurs at a depth of about 25 
g cm~ 2 . Measurements and detailed simulations of the cosmic-ray dose rate have 
been carried out by the Martian Radiative Environment Experiment (MARIE) on 
the Mars Odyssey Orbiter (see Saganti et al., 2004), with measured and modeled 
surface values mostly around 0.06 Gy yr _1 . Integrated over a microorganism's 
lifetime, the corresponding doses are smaller than the lethal dose for any terrestrial 
microorganism, but the whole-body human dose would be significant over decades, 
suggesting problems for manned Mars missions. The dose due to ionizing radia- 
tion is also of interest concerning natural transfer of microorganisms between solar 
system bodies, in particular Mars and Earth (Mileikowsky et al., 2000). 

X-rays from the most energetic solar flares are potentially important for Mar- 
tian habitability, but have not yet been treated. The steady solar coronal X-ray, 
EUV, and FUV emission is currently fairly small (see Giidel et al., 2003). In con- 
trast, during a solar flare, the flux of ionizing photons can increase by orders of 
magnitude, and large dose rates are possible. Unlike on Earth, the Martian atmo- 
sphere is currently thin enough that a significant fraction of this radiation can arrive 
at the surface, where it may lead to genetic damage and mutations in terrestrial-like 
organisms. The fraction of energy absorbed in the atmosphere or reprocessed to 
UV (see Smith et al., 2004a,b) is minimal. Moreover, the X-ray flux from the flares 
of interest here are orders of magnitude larger than the Martian keV X-ray fluxes 
observed by the Chandra X-ray Observatory (Dennerl, 2002) that are probably gen- 
erated in the Martian atmosphere by solar X-rays (see Cravens and Maurelis, 2001) 
and the solar wind (see Gunell et al., 2004) 

Evolution of enhanced resistance to the X-ray flares discussed in the present 
work does not seem possible for anything analogous to terrestrial evolution. The 
flare durations are very short (about an hour), and the times between exposures 
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are orders of magnitude larger than organism lifetimes, so no adaptation could 
occur. Thus highly intermittent flare radiation can be lethal in a way that gradually 
increasing UV (or nearly constant cosmic rays) cannot. 

In this paper, we focus on the following question: At the Martian surface, 
for a strong solar X-ray flare of given total energy and spectrum, how does the 
biological dose compare to the lethal doses of various terrestrial organisms? We 
then extrapolate the recurrence frequency-energy release relation for solar flares 
to estimate the mean time between such lethal events for a variety of organisms. 
Finally, we contrast X-ray dose rates with that of cosmic rays and discuss shielding 
considerations. 

2 Solar Flare Properties 

2.1 Spectra 

X-rays and y-rays with energies between 10 _1 and 10 6 keV 1 comprise a significant 
fraction of solar flare photon emission (e.g. Haisch et al., 1991; Hudson, 1991; 
Kanbach et al., 1993; Ryan, 2000). Hard X-ray (hereafter HXR) spectra are often 
taken to peak around 10-25 keV (Crosby et al., 1993; but see Krucker and Lin, 
2002), but this is very uncertain because it is near the sensitivity cutoff of most pre- 
RHESSI instruments (see Battaglia et al., 2005). Flare X-ray spectra are often fit by 
a single power law (Crosby et al., 1993; Bromund et al., 1995; Veronig et al., 2002) 
or two piecewise power laws (Krucker and Lin, 2002). The range of the estimated 
power-law spectral index (log-log slope) p is large for hard X-ray flares (see Lee 
et al. 1993, Bromund et al. 1995; also Qiu et al. 2004 for HXR microflares): most 
flares have a p between 2.5 and 6, with a median around 4 and little correlation 
with total X-ray output. Using results from RHESSI, Battaglia et al. (2005) find 
that the non-thermal (hard X-ray) emission spectral index is correlated with total 
[ 1 A= 12.4 keV 
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energy release, such that the larger flares have flatter (harder) spectra. These were 
relatively low-energy flares, but if a similar relation holds for the most powerful 
flares, then the large solar X-ray flares described here may deliver even larger doses 
than we estimate (see £|4.2t . Given the uncertainty described, we take the slope of 
the incident photon spectrum as a parameter (see £j3jl, constrained by the observed 
distribution of p, and assume no correlation between spectral index and total energy 
output. As mentioned above, this likely gives a lower limit on the dose estimates. 

2.2 Durations 

Solar flares have a complex temporal structure. The distribution of flare durations 
has been estimated using several datasets and characteristically fits a decreasing 
power law (e.g. Lee et al., 1993). Most HXR flares have short durations (< 20 
sec), with a loose correlation between duration and total energy release. We are 
only concerned with the rare, highest-fluence flares, however, with durations found 
by Crosby et al. (1993) for HXR and Veronig et al. (2002) for SXR flares of 
about 15 min, with the more energetic flares tending to last longer. The X-class 
data of Veronig et al. have a median duration three times longer than the B-class 
data (30 min and 10 min, respectively). Just as with spectral slope, the flare-to- 
flare variation in duration is large and difficult to characterize quantitatively. For 
example, some flares emit much of their energy as gamma rays over several hours 
(see Ryan, 2000). 

For our purposes, though, the duration is unimportant; the acute lethal dose is 
almost always independent of dose rate (Sparrow et al., 1967), and thus depends 
on fluence, not flux. Note that almost all flares have durations much smaller than a 
Martian day, suggesting that mutation and lethality could occur only on one hemi- 
sphere at a time. On the other hand, flare events clustered over many days are 
sometimes observed (e.g., the Halloween flares of Oct-Nov 2003; Woods et al., 
2004) that could progressively irradiate the whole planet. 
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2.3 Energy Release 

The total photon energy release in flares is difficult to estimate and varies by at least 
10 s from flare to flare, but a large number of studies using EUV, soft X-ray (SXR), 
or hard X-ray (HXR) satellite events roughly agree that the differential distribution 
of flare energy releases dN/dW is a power law with index about —1.6 to —1.8 over 
at least six orders of magnitude in total energy release W (see Hudson, 1991; Lee et 
al., 1993; Crosby et al., 1993; Bromund et al., 1995; Aschwanden et al., 2000; Lin 
et al., 2001; Giidel et al., 2003; Qiu et al., 2004; and references therein), although a 
somewhat steeper index (—2.0) has been inferred from a very large sample of SXR 
flares (Veronig et al., 2002). 

Eleven large X-class flares occurred during the extraordinary solar outbursts 
between 18 October 2003 and 5 November 2003, with SXR releases in the GOES 
1-8 A(~ 2-10 keV) band peaking at about 2 x 10 31 erg for an effective duration 
of 30 min. Observations using the SORCE instrument's Total Irradiance Monitor 
yielded a total flare energy at all wavelengths for the somewhat weaker 28 Octo- 
ber flare of 4.6 x 10 32 erg (Woods et al., 2004). Radiation and charged particles 
from these flares compressed the Earth's Van Allen belt to within 20,000 km of 
the surface (Baker et al., 2004), temporarily damaged the orbiting Mars Odyssey 
communication instruments, and significantly reduced stratospheric ozone levels 
(Randall et al, 2005). 

Based on astronomical observations, we believe it reasonable to infer that much 
more energetic solar flares have occurred in the past: 

1. There is little doubt that solar-like stars can produce very energetic, and 
frequent, X-ray flares. Data presented by Audard et al. (2000) indicate SXR flares 
(0.01-10 keV) of energy greater than 3-5 x 10 34 erg for the young, solar-like stars 
47 Cas and EK Dra (age ~ 100 Myr) and 2 x 10 33 erg for an older, solar-like star 
K Cet (age ~ 1 Gyr), with frequencies of about one per ten days at these ener- 
gies. Considering the short sampling period and the form of the derived frequency- 
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energy release relation, more energetic flares seem likely, though less frequent. 

2. Schaefer et al. (2000) have identified nine "superflares" with energy outputs 
of 10 33 to 10 38 erg on otherwise normal F8-G8 main sequence stars. Two of these 
stars were solar-like and produced X-ray flares with energies about 100-1000 times 
larger than the largest observed solar X-ray flare. These flares cannot be attributed 
to binaries, rapid rotation, or youth, and may in fact be common in solar-type 
stars, although cannabalization of giant planets has been suggested (Rubenstein 
and Schaefer, 2000). 

3. Stothers (1980) interpreted NO^ abundance spikes in several Antarctic 
ice cores as due to flares with energies around 10 32 -10 33 erg with a recurrence 
timescale of ~ 10 2 yr, consistent with the extrapolation of the frequency-energy 
relation we adopt below in Eq. |6] 

4. On lower mass dMe main-sequence stars, flares are more frequent, with 
photon releases exceeding 10 34 -10 35 erg observed, often serendipitously (Hawley 
and Pettersen, 1991; Pagano et al., 1997; Liebert et al., 1999; Favata et al., 2000; 
Christian et al., 2003). These energies are enormous, considering that the bolomet- 
ric luminosities of these stars are much smaller than that of the sun. 

Although flares of such large energy release have not been observed in the sun, 
these data for other stars strongly suggest that the maximum observed is limited 
only by duration of observation. We have found no compelling physical argument 
for an upper limit to flare energy releases below the levels discussed in this pa- 
per, but the slope of the empirical frequency-energy release relation suggests that 
some upper limit must exist to avoid energy divergence (Hudson, 1991). Lim- 
its on energetic particle fluxes (e.g., Lingenfelter and Hudson, 1980) are due to 
self-confinement effects and cannot be applied to photons. Given the above con- 
siderations, we expect the the maximum flare release to be at least 10 34 erg and 
assume this to be true in what follows. 
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3 Radiative Transfer Method 

We calculate in detail the transfer of the incident ionizing radiation through the 
Martian atmosphere using a previously developed Monte Carlo code (Smith et al., 
2004b) that accurately treats Compton scattering and X-ray photoabsorption. We 
use the empirical approximation of Setlow and Pollard (1962) for the photoabsorp- 
tion cross section for energies greater than the K edge: 



where E is the energy in units of 511 keV and Z is the atomic number of the 
target atom. We have compared this with the NIST database of cross sections and 
find good agreement. Compton scattering cross sections are calculated using the 
full Klein-Nishina formula. We assume that the atmospheric gas density falls off 
exponentially with height, with a scale height of 1 1 km and a total column density 
of 16 or 22 g cm -2 , corresponding to rough limits set by the freezing and thawing 
of the polar ice caps. This model works well because the high-energy radiative 
processes involved here are relatively insensitive to the exact vertical distribution 
and composition of gas, instead depending mainly on the total column density and 
mean molecular weight. 

Solar flare photon number spectra are assumed to be distributed as E~ p , with 
2 < p < 6 (see Sj2ji. For calculational purposes, the flare spectrum is assumed to 
extend from 10 keV to 511 keV. The 10-keV lower limit is taken because photons 
below this energy are nearly completely absorbed by the atmosphere and unimpor- 
tant to the surficial dose, while the upper limit is set high enough that a negligible 
number of incident photons are at higher energies, even for the shallowest spectra 
(lowest p). Given a lower limit E m { n and a specified normalization Af tot , the differ- 
ential number spectrum can be written (neglecting any upper cutoff for simplicity) 




(1) 




(2) 
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and the average energy in the spectrum is obtained by integrating between the 
minimum and maximum cutoff energies: 



This is the spectrum that is propagated through the atmosphere using the Monte 
Carlo code. Note that henceforth the total energy release figures refer only to the 
energy between 10 and 511 keV. 

After transport through the atmosphere, the photon fluence is converted to a 
biological dose by assuming that the surficial radiation is absorbed by pure water. 
Much of the damage by ionizing radiation is thought to be "indirect," involving 
chemical reactions initiated by energy deposited in the bulk cell water or first hy- 
dration layer rather than "direct" ionization of DNA (von Sonntag, 1987; see Ward, 
1999, and references therein), although this terminology is now recognized as an 
oversimplification (see Fielden and O'Neill, 1991). We then estimate the dose 
by integrating the surficial energy spectrum E dN/dE over the energy-dependent 
opacity of water K W (E): 



The water opacity in our code includes both Compton scattering and photoabsorp- 
tion. It should also be noted that the biological quality factor and relative biological 
effectiveness (RBE) are both near unity for X-rays and y-rays, so no adjustment is 
required to compare to empirical lethal doses, as in cosmic-ray studies. 

Most of the uncertainty in our calculations is due to the fact that the total X- 
ray opacity of the Martian atmosphere depends sensitively on the uncertain energy 
spectrum of the flare photons, since the photoabsorption optical depth varies ap- 
proximately as exp(— E~ 3 ). Given a particular flare spectrum, our results are ac- 
curate, but the spectral slopes vary greatly from flare to flare, even given the same 
total energy release (see Thus we use the slope, p, as a free parameter that 
varies between 2 and 6 and present the results for a range in p. Later we integrate 
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over an approximate representation of the observed frequency distribution of p for 
hard X-ray solar flares to derive appropriate average values for waiting times as a 
function of dose. 

4 Results 

4.1 Calculated Surficial X-Ray Spectra 

The Martian atmosphere, though thin by terrestrial standards, significantly atten- 
uates the incident flare radiation, but the high fluences associated with the most 
massive flares lead to extreme irradiation of the surface nonetheless. We show in 
Fig. [^the surficial fluence spectra from a 10 35 erg flare for two atmospheric col- 
umn densities (E = 16 and £ = 22 g cm~ 2 ) and two spectral indices (p = 4 and 
p = 2.5), together with the above-atmosphere, incident spectra. The vertical axis 
scales linearly with flare energy release (we have used 10 35 erg for this example). 
In general, for all column densities, spectra with lower mean energies (higher p) 
are attenuated more because of the rapidly rising (E~ 3 ) photoabsorption opacity. 
This effect, combined with the falloff of the incident spectrum at high energies, 
together creates a peak in the surficial spectrum at moderate energies (30-50 keV 
here). Additionally, increasing the atmospheric column increases the attenuation 
at all energies, as one would expect. 

4.2 Estimated Biological Doses 

Ionizing radiation alters DNA through direct ionization and indirect chemical inter- 
actions involving diffusion of radical products of ionization and dissociation in the 
cytoplasm. Resulting biochemical effects include single- and double-strand breaks, 
base damage or abasic sites, multiply damaged sites, oxidized base clusters, and 
cross-linking within DNA or with proteins (see Becker and Sevilla, 1993; Ward, 
1999). An amazing suite of repair processes and associated enzymes repair the 
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Figure 1: Surficial, reprocessed, X-ray flare spectra, assuming a parameterized, 
power-law model spectrum (see Two incident spectra are shown, with power- 
law spectral indices of p = 2.5 and p = 4. For each of these spectra, associated 
surficial spectra are shown for two different column densities: 16 and 22 g cm~ 2 . 
Attenuation due to photoabsorption is most severe at low energies, leading to a 
sharp decline below about 20 keV in the flux per unit energy. Also, the original 
incident spectrum falls off at high energy, so the combined effect leads to a peak at 
intermediate energies (30-50 keV here). 
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damage (Friedberg, 2003, and references therein), but generally for each organism, 
above some critical amount of absorbed energy, fecundity greatly diminishes. Bi- 
ological damage is quantified in terms of the amount of energy absorbed, or dose, 
typically in units of rads or Grays (1 rad = 100 erg g _1 absorbed; 1 Gray = 100 
rad). The lethal dose is often defined as the dose that kills a particular fraction of 
a laboratory population. The dose for reducing the population by l/e is denoted 
D37, while the dose for reduction of the population by 90% is D90. 

The calculated X-ray doses in our situation are shown in Fig. |2] as a function 
of power-law index for two column densities: £ = 16 g cm~ 2 and £ = 22 g cm -2 . 
The calculation assumes a solar flare with a total X-ray release of 10 35 erg (see ^2]i. 
Again, for flares of different energy release, the dose simply scales proportionally. 
The doses for the shallowest spectra (2 < p < 3) are about two orders of magni- 
tude higher than those for the steepest spectra (5 < p < 6), showing the tremendous 
sensitivity of the dose on the form of the incident spectrum. For comparison to the 
power-law spectra, we also show monoenergetic spectra with energies equal to the 
average energy in the power-law spectra for each spectral index. The dramatic de- 
cline in the dose with energy for the monoenergetic cases illustrates the importance 
of the high-energy tail in the power-law spectra and the inaccuracy of doing this 
calculation using a monoenergetic approximation. 

The results shown in Fig. |2]can be confirmed by considering a 10 35 -erg flare 
incident on the top of the atmosphere and subject to only absorption due to Comp- 
ton scattering and photoabsorption. Strictly speaking, the Compton process is a 
scattering process, not absorption, but we will neglect photons that backscatter in 
the atmosphere and eventually reach the surface in order to get this estimate. This 
enables us to write an integral that approximates the total dose D in terms of the 
atmospheric optical depth x a {E) and the water opacity K W (E): 



where N(E) is the incident number spectrum at the top of the atmosphere and E, 




(5) 
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Figure 2: Surficial X-ray dose for a flare with 10 35 erg total energy release versus 
spectral index. (Doses scale linearly with energy release.) Two column densities 
are used: 16 and 22 g cm~ 2 . Power-law flare spectra are denoted "power" in the 
legend. Monoenergetic incoming photons with energies equal to the mean incident 
flare energy are shown as "mono." This case demonstrates the inaccuracy of ig- 
noring the full radiative transfer, since the curve decreases much more rapidly than 
the power-law cases. The final pair of curves (denoted "abs est") uses a power-law 
incident spectrum with the radiative transfer done using the analytic, absorption 
approximation discussed in ^4.21 The general trend among all curves is that the 
dose decreases with increasing spectral index. This is because (as can be seen on 
the upper horizontal axis) the average energy of the flare spectrum decreases with 
increasing spectral index. The atmosphere is vastly more opaque to lower-energy 
X-rays, which dramatically affects the surficial dose. The right-hand vertical axis 
shows tick marks at the mean-log lethal dose of a few representative terrestrial 
organisms (see Table 1). 
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is the minimum energy considered in the Monte Carlo code. This integral can be 
evaluated numerically to confirm the results given in our Monte Carlo calculation; 
we have done this and present it as the final set of curves in Fig. |2] marked "abs 
est." The integral is very close to the full calculation at all energies, being only 
different by at most a factor of two, which suggests that backscattered radiation is 
only marginally important to the surficial doses. For the flattest spectra (p near 2), 
the code gives a higher (and more accurate) dose because it includes backscattered 
radiation. For steeper spectra, the difference between the above approximation and 
the code stems from the discrete number of bins in the surficial spectra in the code 
and the fact that the mean incident spectrum energy is very close to the lower cutoff 
energy, causing discretization errors in the output spectra when binning photon 
energies close to the cutoff. The integral above, being continuous, suffers no such 
errors, so is probably more accurate in this regime, where absorption dominates 
anyway. 

The horizontal ticks along the right side of Fig. |2]show typical mean-log acute 
lethal doses for a variety of organisms. The most comprehensive compilation we 
know of, Sparrow et al. (1967), gives D90 for X-ray and y-ray exposure of 79 
organisms. We have expanded and supplemented this with more recent work, es- 
pecially for microorganisms, including the review of virus lethal doses by Rohwer 
(1984), lethal doses in bacterial spores summarized by Russell (1982), studies of 
food pathogenic bacteria by Dion et al. (1994) and Farkas (1998), studies of Ar- 
chaean hyperthermophiles by Kopylov et al. (1993) and DiRuggiero et al. (1997), 
mammalian studies presented in Fielden and O'Neill (1991), studies of mutations 
in human survivors of Hiroshima and Nagasaki (see Turner, 1995, p. 398), the 
Chernobyl accident and Kazakhstan nuclear weapon tests (Dubrova, 2003), and 
inhabitants of the Kerala (India) radioactive hotspot (Forster et al., 2002). We find 
these results to be consistent with each other in cases we have checked. 
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Table 1 : Acute lethal dose ranges for representative organisms 

Organism Acute Lethal Dose (Gy) 

mammals 1-10 

higher plants 5-700 

insects 1-2000 

bacteria, protozoa, algae 30-15,000 

viruses 150-20,000 

bacterial spores 1,000-4,000 



Because of the large range of lethal doses found in nature we have chosen to 
represent various organisms in Fig. EJby the average of the base- 10 logarithm of 
the lower and upper limits found for that class. This is an appropriate choice be- 
cause the available data from Sparrow et al. (1967) within each taxonomic class are 
distributed roughly uniformly in log lethal dose. Thus for mammals, with a range 
of 1-10 Gy, the number indicated on the log scale is 0.5. Table 1 lists the approxi- 
mate lethal dose ranges for these organisms. Organisms with extremely large lethal 
doses are interesting for long-term survival prospects — e.g. adult Drosophilia at 
1000 Gy, some Archaean hyperthermophiles at 2000-6000 Gy, Rubrobacter ra- 
diotolerans at 8000 Gy (Ferriera et al., 1999), Deinococcus radiodurans at 10,000 
Gy (depending on cell phase; see Battista, 1997), and the virus-like proteinaceous 
infectious particle (prion) associated with scrapie at 20,000 Gy (Rohwer, 1984). 

Even very low doses can be important. Doses for mutation and clustered 
DNA damage can be much smaller than the lethal dose. For example, in humans, 
the whole-body acute lethal dose is about 1-2 Gy, but mutations (Sparrow et al., 
1972; Sankaranarayanan, 1982; Forster et al., 2002), chromosomal abnormalities 
in blood lymphocytes (Violot et al., 2005), and clustered DNA damage (Suther- 
land et al., 2000) occur at only around 0.01-0.2 Gy, and germ-line mutation dou- 
bling is believed to occur around 0.5 Gy (UNSCEAR 2001; Sankaranarayan and 
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Chakraborty, 2000; but see Forster et al., 2002). So doses even smaller than the 
"humans" level in Fig. |2]should be considered hazardous for human exploration of 
Mars. We indicate this level at 0.1 Gy with the label "human risk" and emphasize 
that this value is illustrative, not definitive. 

We have ignored the dependence of lethal dose on dose rate. Sparrow et al. 
(1967) claim no evidence for a significant dependence in the literature they com- 
piled, but it is undoubtedly true that at some level a dose rate dependence occurs, as 
has been observed in mammalian and other organisms. And since the duration of 
solar flares is random, it is difficult to convert received fluences into fluxes, which 
are needed for dose rate calculations, hence we have ignored this and focused en- 
tirely on time-integrated quantities in this work. 

Finally, recall the possible correlation between hard X-ray flux and spectral 
index found in the RHESSI data by Battaglia et al. (2005) mentioned in £12.11 If 
flares with a larger total energy release do indeed typically have flatter spectra, then 
they would have higher mean photon energies. The Martian atmosphere is more 
transparent to X-rays of higher energies, so this effect would boost the surficial 
doses received from larger flares. So in this sense our results are an upper limit to 
the flare releases required to produce given doses. 

4.3 Mean Times Between Lethal Events 

The calculated surficial doses are very sensitive to the flare spectral index p, which 
varies greatly and unsystematically with other flare parameters, so a better way 
to determine the lethality risk is to integrate the frequency of events delivering a 
given dose for each p over the probability distribution of p values. We take this 
probability distribution to be the uniform distribution, which is roughly consistent 
with results shown in Lee et al. (1993), Bromund et al. (1995), Veronig et al. 
(2002), and Qiu et al. (2004) and other references given there. This procedure 
yields the total frequency of lethal events for an organism at the Martian surface 
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due to all flares. 

To estimate flare timescales as a function of release W, we normalize to one 
10 32 -erg solar flare per decade. Extrapolating the average HXR frequency-energy 
statistic dN/dW (see Sj2ji gives the average time between events of energy release 
W32 (in units of 10 32 erg) as T(W) ~ 10W 3 2 yr, with q between about 1.6 and 1.8, 
and the mean time between events at least as large as W32 as 

T(>W 32 ) = \0(q-l)W 3 " 2 l yr. (6) 

This agrees with the estimate of Hudson (1991) and is broadly consistent with the 
dozen or so 10 31 -10 32 erg events that have been observed since GOES soft X-ray 
monitoring began in 1976. Upper limits on proton fluences in lunar rocks, tree ring 
records of 14 C, and the requirement of finite energy suggest that the frequency- 
energy release relation dN/dW must steepen above some energy (Reedy et al., 
1983; Hudson, 1991; Aschwanden, 1999), but the value of that energy is unknown. 

Figure |3] graphically shows the mean time between lethal events for the same 
set of organisms given in Fig. |2]as a function of surficial dose, for an atmospheric 
column density Z of 16 g cm~ 2 . To produce this plot, we use the mean of the 
empirical values of frequency-energy release slope, so q = 1 .7 in Eq. |6] The right- 
hand side vertical axis shows the corresponding flare energy releases W required 
for lethality. Again, this result is integrated over the spectral index distribution. 
Additionally, for a column density of 22 g cm~ 2 , the attenuation is roughly a factor 
of 2-3 higher (as can be see in Fig. 12), so the required lethal dose would be roughly 
a factor of 2-3 higher, implying a mean waiting time that is 2 0J -3 0J = 1.6-2.2 
times longer. 

The mean time to lethality is strongly weighted toward the flattest (p ~ 2) 
flares. The atmosphere is more transparent to the flatter flares (because they have a 
higher average energy), so a lethal fluence at the surface requires less total energy 
release at the sun. Since flares with lower total energy release are much more 
common, chances are that the flare that finally delivers the lethal dose will have a 
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Figure 3: Mean time between flare events that deliver a given dose. The dose deliv- 
ered depends on the spectral index, such that flares with larger spectral index must 
have a larger total energy to deliver a given dose (see Fig. |2j. We have integrated 
over the distribution of spectral indices, and accounted for the dose dependence on 
this index, in order to find a mean time between events that deliver a given dose 
without specifying the spectral index (see £)4.3t . The upper horizontal axis shows 
the lethal doses of a variety of organisms for reference (see Table 1). As expected, 
the time between events increases as the dose necessary for lethality increases. 
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smaller spectral index. We show this explicitly below. 

To derive the equation for the mean time between flares of a given dose, as 
shown in Fig. |3j we first convert the waiting time distribution given in Eq. |6]back 
to the empirical flare energy release-mean frequency relation: 

v( W ) = k w W- 2 q yr-\ (7) 

where once again we use q = 1.7 and we estimate kw ~ 0.1. The flares are roughly 
distributed uniformly by spectral index p, and since p ranges from 2 to 6, we 
can write the probability distribution function of p as f(p) = 1/4. The frequency 
distribution for the flares as a function of energy W and spectral index p is then 

v(W,p) = v(W)f(p), (8) 

where we have assumed that W and p are uncorrected. 

To predict how often solar X-ray flares can deliver a particular lethal dose, we 
must transform this distribution to a function of dose D. This requires knowledge of 
the relation between total flare energy release and surficial dose for each spectral 
index p. In our case, this relation comes from our radiation transport code. An 
empirical fit to our calculations for a column density of 16 g cm~ 2 shows that the 
dose in Gray as a function of flare energy is 

D(W,p) = — e-P/P° Gy, (9) 
k D 

with ko = 0.23 and po = 0.53. We can write the energy required to deliver a given 
dose by inverting this equation. Then, using the standard technique for transform- 
ing probability distribution functions, we have 

d(W,p) 



v(D,p) = v[W(D,p),p 



v[W(D,p)\ 



d(D,p) 
dW 



(10) 



dD 



f(p) (11) 



k -^l D -^xp[(l- q )p/p ] yr- 1 . (12) 
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Next, we would like to know how many flares of all spectral indices occur each 
year, but now the dose and spectral index are not independent, so we must integrate, 
instead of just scaling by 1/4: 

r 6 



v(D) = J v(D,p)dp (13) 
= ^_ g 1 f7 ) ° £>- g {exp[2(l -q)/po] -exp[6(l -g)/po]} y^ 1 . (14) 



All flares that deliver doses above the lethal dose are most definitely lethal, so 
we must integrate the flare dose-mean frequency distribution above the lethal dose 
(we use D90 here) to find the total number of flares that occur on average per year 
that deliver at least the lethal dose: 



v(>D 90 ) = / v(D)dD (15) 

= k ^t0 { ex P[2(! " l)/Po] ~ «p[6(l - q)/po}} yr-^6) 

Finally, to get the mean time between lethal events, simply take the inverse of this 
frequency: T(> D90) = 1 /v(> D90). For the case discussed in this work, kyy ~ 0.1 
and q = 1 .7. Thus the typical time between lethal flares as a function of lethal dose 
in Gray is 

T(>D 90 )~ 200 D° 7 yr. (17) 

This time is much shorter than an estimate using just a particular p value yields. 
For example, using p = 4 instead of integrating over the distribution of p would 
give about 3000 yr for the coefficient above instead of 200 yr. Properly treating 
the spectral indices as a distribution in the calculation is important, rather than 
just assuming that p = 4 is the representative index and making simple estimates 
based on that. In fact, the true waiting times are probably a bit shorter than our 
estimate, since the distribution of spectral indices isn't exactly flat, but rather is 
peaked slightly at lower indices. In other words, the average flare spectrum may be 
somewhat harder than we take here. 
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Figure |3] indicates that a solar X-ray flare of 10 38 erg would be needed to kill 
any organism as resistant as D. rad, bacterial spores, and the most resistant viruses 
(see Table 1) on the surface of Mars. Such flares would occur at least every ~ 
200 kyr if our extrapolation of the solar frequency-energy release scaling could be 
extended to these extreme energies. Since such energies have not been observed 
in the sun or other (non-binary) stars, we conclude that such extremely radiation- 
resistant microbial life on Mars may not be affected by solar flares. Bacteria and 
protists with the adopted mean log lethal dose require slightly less extreme flare 
energies and have corresponding times to lethality of about 20 kyr. Reptiles, in- 
sects, and plants require more reasonable flare energies, especially considering that 
many have lethal doses much less than the mean log (see Table 1 and Sparrow et 
al., 1967). Flares of only 10 34 — 10 35 erg will result in mutations and lethality for 
mammals, reptiles, insects, and plants, and may occur every few hundred years. 
Similarly, a 10 34 -erg flare, which can be lethal to humans, might occur every 200 
yr on average. We conclude that no organisms other than the most resistant bacteria 
and viruses could have survived on the surface of Mars during the past ~ 10 5 yr, 
which is almost certainly much smaller than the time during which Mars has had 
its current low-column-density atmosphere. 

Additionally, a correlation between hard X-ray flux and total flare energy re- 
lease in which larger flares have harder (flatter) spectra (e.g., Battaglia et al., 2005) 
would decrease the time between lethal events. Larger flares would be more pen- 
etrating and would deliver higher doses than we have estimated, which would de- 
crease the energy release needed for a given dose. Since less energetic flares are 
more common, lethality would occur more often. 
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4.4 Shielding and Dose Rate Considerations 
4.4.1 Soil and C0 2 Ice Shielding 

So far, we have ignored possible protection strategies, such as radiation shielding 
sources. As an illustrative example, we take the Mars Pathfinder rover soil sample 
data of Wanke et al. (2001) and calculate the opacity of this material to incom- 
ing flare spectra. Based on the composition of the soils in Table I of Wanke et 
al. (2001), we find that the samples have an effective atomic number of 15.2 (the 
Z that, when used in Eq. ^ yields the mean cross section of the soil obtained 
by element-by-element addition). The combined photoabsorption and Compton 
scattering opacity of this material can be calculated. We use the analytical approx- 
imation in Eq. |5]to calculate the transfer of the surficial radiation through the soil. 
Using an atmospheric column density of 16 g cm" 2 , we find l/e (i.e., optical depth 
unity) shielding columns for X-ray flares of spectral index p = 2 and p = 6 to be 
7.7 g cm" 2 and 3.6 g cm -2 , respectively. For an atmospheric column density of 
22 g cm" 2 , the same calculation gives only a slightly different result: 8.1 g cm -2 
and 3.9 g cm -2 , respectively. With a typical density of terrestrial andesites of 2.7 
g cm" 3 , these soil columns suggest that around a few cm of regolith is enough to 
provide l/e attenuation and that tens of cm will reduce the dose by 99% or more. 
For comparison, the possible ice layers found by Boynton et al. (2002) at 40-150 
g cm -2 below the surface should be very well shielded. 

CO2 ice has an effective Z of 7.4, roughly half that of the soil, so l/e shielding 
requires a slightly larger ice column. For an atmospheric column density of 16 g 
cm" 2 , we find l/e shielding columns of CO2 ice for spectral indices p = 2 and 
p = 6 to be 8.2 g cm" 2 and 4.6 g cm" 2 , respectively. For a column density of 22 
g cm" 2 , we get 8.6 g cm" 2 and 4.9 g cm" 2 , respectively. Roughly three times the 
thickness of ice would be needed for the same l/e shielding as corresponding soil 
columns due to the different densities. The difference is smaller than expected from 
the Z 3 dependence of the photoabsorption cross section because a large fraction of 
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photons in the surficial spectrum, even for steep flare spectra, has energies high 
enough that the radiative transfer is dominated by Compton scattering, a process 
insensitive to elemental composition. 

Altitude can change the atmospheric column density and have a measurable ef- 
fect on the dose. The column density above Hellas Basin, the lowest topographical 
point on Mars, is roughly 8 to 14 g cm~ 2 larger (depending on season) than the 
average and consequently reduces the dose to a factor of about 1/7 to 1/3 of the 
dose at the mean surface level. This suggests that organisms in low areas would be 
better protected, but not completely safe. Hellas Basin is an extreme example, and 
more generally we expect topographic altitude variations to give less than a factor 
of two variation in dose for a given flare. 

4.4.2 Manned Mission Shielding 

Astronaut protection is a critical issue for long-term, manned space missions. The 
total mass of a spacesuit depends mostly on the column density of the component 
material, so, given a fixed shielding requirement, the suit's mass will be minimized 
by maximizing the suit's opacity to incoming radiation. Typical next generation 
spacesuit designs provide only up to a few g cm~ 2 of mostly low-Z shielding mate- 
rial (see Ross et al., 1997, for a review). The lightweight and mostly low-Z suits are 
designed to provide efficient protection from incoming energetic particle radiation 
(such as cosmic rays). Charged particles lose more energy per unit column density 
traversed in materials with lower Z because the energy loss is roughly proportional 
to Z/A, where A is the mean atomic mass, and lighter elements maximize Z/A. 
Additionally, low-Z materials produce less secondary radiation (mostly neutrons) 
and lower absorbed doses when incident particles cause nuclear reactions in the 
shielding material (see, e.g., Schimmerling et al., 1996, and Wilson et al., 2001). 

X-ray radiation is most efficiently stopped by high-Z materials. Neglecting 
engineering considerations and particle radiation, the ideal shielding material for 
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X-rays would maximize the quantity Z 3 /A, because the photoabsorption opacity, 
which dominates at these energies, has roughly this dependence on Z and A, as can 
be seen in Eq. [2 Thus, heavier elements are favored for X-ray shielding, since 
atomic mass A increases much slower than Z 3 across the periodic table. Spacecraft 
and surface structures, such as "Marsbases," would probably contain a significant 
amount of high-Z material, and thus could provide adequate flare X-ray shielding 
if surrounding structures provide a shielding column of greater than tens of g cm -2 
of Z > 10 material. For example, polymers, water, and lower-Z metals, such as 
magnesium and aluminum, would be less desirable than titanium or steel alloys 
(especially those high in molybdenum). And even a few g cm~ 2 of high-Z ma- 
terial is inadequate should a large solar X-ray flare occur during a spacewalk or 
"Marswalk." Balancing the two contrasting shielding requirements for particles 
and photons while minimizing the weight and rigidity of spacesuits may prove to 
be difficult. 

All of these considerations also apply to moon missions and space missions 
elsewhere in the solar system, with doses being even higher outside of Mars' pro- 
tective atmosphere and inversely proportional to the square of the distance to the 
sun. 

4.4.3 X-ray Flares and Cosmic Rays: Dose Rate Differences 

The mean times to lethality we give refer only to the surface of Mars on the hemi- 
sphere facing the sun for organisms with short generation times. Lethality from 
a solar flare contrasts with the situation for cosmic-ray exposure; the difference is 
one of acute versus extended weak exposure. Flares are intermittent bursts, and 
lethality will occur during the flare if the flare is sufficiently energetic. The or- 
ganism lifetime does not affect lethality unless it is less than the mean flare dura- 
tion (about 10-60 min for the flares considered here), and even most bacteria have 
longer generation times. For example, the average generation time for many E. coli 
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strains is about an hour. Cosmic rays, however, cause damage over long periods 
of time from constant, relatively small, surficial fluxes. Organisms might be able 
to continuously repair the cosmic-ray damage (as on the present Earth), while the 
high dose rates of flare exposure may overwhelm genomic repair systems. The 
Saganti et al. (2004) and Pavlov et al. (2002) calculations cited in ^give Galactic 
cosmic -ray dose rates of 0.06-0.2 Gy yr _1 . Taking the representative dose rate to 
be 0. 1 Gy yr~ 1 , the time to accumulate a lethal dose from cosmic rays is longer 
than the lifetime of all but the most sensitive mammals and long-lived plants. 

Pavlov et al. (2002) also suggest that living organisms would survive cosmic- 
ray irradiation because of their short lifetimes, and only dormant organisms, such 
as spores, with their inactive genomic repair systems, would be killed by the slowly 
accumulating damage. In contrast, solar flares with very high energies have dose 
rates high enough to kill living organisms (because the dose is so acute) and pos- 
sibly contribute to or even dominate the death of dormant organisms. Dormant 
organisms have radiation tolerances generally much larger than active organisms, 
but they accumulate damage over a much longer time period. In this case the total 
dose would be due to a steady dose rate of cosmic rays plus the accumulated acute 
doses of a large number of individual X-ray flares over the period of dormancy. 

5 Summary and Implications 

Continuing astronomical observational programs are recording larger and larger X- 
ray flares from stars of solar mass and smaller. The most energetic flares observed 
on the sun are almost large enough to cause mutation and lethality in higher organ- 
isms on the surface of Mars, and extrapolation to higher energies implies greater, 
but less frequent, biological damage. 

Using a previously developed Monte Carlo radiative transport code (Smith et 
al., 2004b), we have calculated the surficial X-ray spectra and corresponding bi- 
ological doses due to solar flares with a range of spectral indices. We compared 
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our results to the average lethal doses for a representative range of terrestrial or- 
ganisms. We find that, if the sun has flares much larger in energy output than those 
observed (as supported by observations of other stars), organisms with a range of 
radiological tolerances, from mammals to plants, insects, and reptiles to even high- 
tolerance bacteria, spores, and viruses could be killed on the hemisphere of Mars 
facing the sun during the flare. Though a single flare leads to only partial irradia- 
tion, flares are often clustered in time, so it may be possible to irradiate the whole 
surface over the course of weeks during a period of extreme solar activity. 

Estimates were given for the mean time between lethal events based on ex- 
trapolating the observed recurrence frequency-energy release distribution for solar 
flares. Precise results are difficult because of the random nature of flare occurrence 
and spectral indices, so we provide mean times between lethal flare events inte- 
grated over the spectral index distribution. In addition, the estimated mean times 
between flares of such energies are orders of magnitude less than the age of the so- 
lar system, so may have occurred 10 3 -10 7 times for the period during which Mars 
has had its present thin atmosphere, depending on the maximum possible solar flare 
energy and required lethal dose. 

Finally, we showed that a column thickness of ~ 4-8 g cm~ 2 of Martian soil 
would be sufficient to provide l/e attenuation of the surficial flare energy. The 
conflicting shielding requirements for particle radiation and the X-rays discussed 
here leads to a difficultly designing safe spacesuits for manned Mars or lunar ex- 
ploration, but designing safe "Marsbases" and other ground structures should be 
easier. 
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